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Abstract 


Predicted rates of fire spread using a mathematical model 
were consistently greater butin reasonably close agreement with 
rates observed on test fires in ponderosa pine and Douglas-fir 
slash. Fuel loading, bulk density, particle density, particle 
surface-to-volume ratio, heat content, total plant salt content, 
silica-free salt, fuel moisture, and wind velocity were deter- 
mined as input variables for the test on plots containing slash 
with and without needles, and at two loadings and three depths. 
Fuel discontinuities and inadequate knowledge of fuel moisture 
contents that limit ignition are discussed as the primary reasons 
for the deviations. One-half of the total fuel weight loss in the 
flame front was accounted for by particles less than 1 cm. in 
diameter. A multiple regression for predicting flame length is 
provided. 
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INTRODUCTION 


Development of a mathematical model for predicting rate of fire spread is part of 
an effort at the Northern Forest Fire Laboratory to formulate a system for appraising 
the fire behavior potential of fuel. Desirable features of such a fuel appraisal system 
are quantitativeness, predictability, and flexibility. 


A fuel appraisal system based upon quantitative inputs and outputs is needed for 
objectively guiding fire management decisions. Appraisal of the fire spread potential 
of fuel can help with evaluations of silvicultural treatments, prevention and presuppres- 
sion planning, budgeting of fire control expenditures, fire simulation training, and 
wildfire suppression. 


Earlier and existing systems of fuel appraisal, of which Hornby's (1935) is the 
best known and most widely used, have been useful, but contain weaknesses that limit 
their application for current and future resource management. The major weaknesses 
seem to involve the following: 


1. A fixed level of fire-weather severity is required for making the evaluations. 
"Average bad" conditions are assumed. Usually fuels are considered quite dry and winds 
rather strong. 


2. Seasonal and yearly changes in the amount and arrangement of fuel cannot be 
accommodated. 


3. Prediction of fire behavior using changing fuel and weather conditions is not 
possible. 


4, The systems are qualitative. They require experienced people to assess the 
fuels providing a precise meaning only to those closely familiar with the fuels. 


A quantitative basis for fuel appraisal, such as provided by a mathematical model, 
facilitates uniform interpretation of appraisals and makes predictability and flexi- 
bility possible. Predictability permits assessment of anticipated land management 
practices upon the vegetation as fuel. For example, the fuel conditions created by 
harvesting and thinning could be predicted and appraised before cutting, and plans 
formulated for proper treatment of the expected hazard. Flexibility means that the 
application of fuel appraisals could be as broad or specific as desired, depending on 
the fuel characteristics used as modeling input and their representativeness. 


A system of fuel appraisal consists of two phases: (1) Measurement and descriptia 
of fuel properties required as input to the system; and (2) rating the fire behavior 
potential indicated by the fuel properties. This paper concerns the second phase; 
specifically, it reports on tests to determine the accuracy of a mathematical rate-of- 
fire-spread model designed by Rothermel (1972). 


THE RATE-OF-SPREAD MODEL 


Rothermel's rate-of-spread model is summarized by the following equation: 
EEG $4 +o)') 
Rate of spread, R oy oA A etl ee (1) 
PVs 
where: 


I, = reaction intensity or the rate of energy release in the flame front of a fire 
B.t.u./ft.2-min. 


— = propagating flux ratio--the proportion of I, required for propagation of a 
no-wind fire, dimensionless 


ee = multiplying factor for the effect of wind, dimensionless 
ie = multiplying factor for the effect of slope, dimensionless 
= bulk density of fuel, 1b./ft.3 


e = effective heating number which expresses the proportion of fuel that is 
heated to ignition temperature dimensionless 


Qe = heat of preignition, the heat required to bring a unit weight to ignition, 
Bet.) 10. 


These terms are computed from a set of rather complex equations. This set, which 
is summarized in equation (1), constitutes the mathematical model. The basic fuel and 
environmental input variables for the model are shown in table l. 


Table 1.--Fuel and envtronmental inputs that determine values for the terms in 
the rate-of -ftre-spread equation 


Input variables Terms in the spread equation 
ee) ae eS" a See en ome (ot er Gig Qig 
Loading, 1b./ft.2 X X X X X 
Heat content, 
Bat eu (bi X 
Fuel particle density, 
Ibe /ft. 2 K X X X 


Fuel particle surface- 
to-volume ratio, 


fteo/ eta X x X X 
Depth: of fuel’, ft’. X X X X X 
Fuel moisture content, 

fraction of dry wt. X X 
Total salt content, 

fraction of dry wt. X 


Silica-free salt 
content, fraction 
of dry wt. X 
Wind velocity, 
ft./min. X 
Angle of slope X 


The model was developed from data gathered using controlled laboratory combustion 
facilities, fuel beds comprised of uniformly distributed dead fuel particles of a single 
size, and quasi-steady-state fire. Conditions were varied to show the effect of fuel 
depth, loading, particle size, windspeed, and slope. The effect of fuel moisture and 
mineral content was taken from other work. 


Predictions of fire spread and flame front intensity using the model should apply 
most accurately to surface fires in uniform distributions of such fuels as needle lit- 
ter, grass, logging slash, and perhaps crowning fires in brush. Less accurate predic- 
tions would be expected for patchy distributions of fuel. 


The model pertains to fire spread by a flaming fire front. Fire propagation from 
spotting, fire whirls, and blowup conditions is not a part of the model. The model 
was designed for fuel appraisal, fire-danger rating, and presuppression planning where 
relative fuel and fire behavior evaluations are useful. Accurate fire behavior predic- 
tions for specific fires will require further refinement of the model and understanding 
of fuells. 


Application of the model to heterogeneous fuels--fuels containing more than one 
size of fuel--requires recognition of the contribution to fire spread made by each kind 
and size class of particle. This is accomplished by entering the input variables as 
mean values weighted by surface area of each kind and size class of fuel particle. A 
complete explanation of the model and its development is given by Rothermel (1972). 


PROCEDURES 


Observed and predicted rates of fire spread were compared using test fires on an 
old bed of the Clark Fork River east of Missoula, Montana (fig. 1). 


Plots 8 by 30 feet were prepared using logging slash hauled by truck as the 
fuel. All other vegetation was cleared from the plots. Slash was used as a test fuel 
because it (a) often is a hazard and silvicultural problem of considerable concern to 
land managers, (b) cannot be readily tested in a laboratory, and (c) provides a test 
of the spread model using heterogeneous fuel for which only theory exists. 


Slash of two species, Douglas-fir (Pseudotsuga menztestt var. glauca (Beissn.) 
Franco) without needles and ponderosa pine (Pinus ponderosa Laws.) with needles, was. 
used at three depths for each of two loadings. Each combination of species loading 
and depth was replicated twice, making a total of 24 plots. Two extra plots were con- 
structed for training but were not needed for this purpose so that 13 test fires were 
burned for each species. The loadings and depths are within a range that is often 
encountered in slash. 


The plots were on a zero slope and were burned under low wind velocities. A test 
of the spread model at high wind velocities was not possible. Fairly low fuel moisture 
contents were sought and obtained on all fires. 


The plots were loaded uniformly with weighed amounts of slash providing conditions 
of approximately 15 and 30 tons per acre. Immediately after weighing and loading, the 
plot moisture content was sampled by ovendrying needles and branches from several 
diameter size classes. The moisture estimates (with standard errors within 4.9 and 8.6 
percent of mean values for Douglas-fir and ponderosa pine, respectively) permitted pre- 
cise estimation of the actual ovendry weight of fuel. The slash on some plots was 
buoyed with wire to provide desired fuel depths. 
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Figure 1.--Plots 8 by 30 feet ie, 
were constructed on an old 
bed of the Clark Fork River. 
Wtre strung between steel 
fenceposts supported the 
slash at specified depths. 


Measurement of Spread Model Inputs 


Loading by fuel particle stze class.--The amounts of fuel in 0- to 1-, 1- to 3-, and 
3- to 5-cm. diameter size classes and the needle class for ponderosa pine were calculated 
from sample determination of the proportions of total fuel in each size class (fig. 2). 


Three branches per plot were randomly collected, dissected into the size classes, 
and the proportions of total branch weight were determined for each size class. These 
proportions were used to calculate fuel weights by size class for ponderosa pine 
because the proportions appeared to be the same on all plots. The amount of fuel in 
each size was actually sampled on each Douglas-fir plot using the planar intersect 
technique (Brown 1971) because the amount of O- to l-cm. size class material appeared 
to vary considerably from plot to plot. 


Parttele surface-area-to-volume ratios.--For the 0- to l-cm. Douglas-fir, an average 
particle surface-area-to-volume ratio (c) was determined by taking random diameter 
measurements. For a cylinder, o is 4 divided by the branch diameter. Diameters for 
the ponderosa pine and the other Douglas-fir size classes were taken at 10-cm. intervals 
along the three sample branches per plot. Standard errors of the estimate averaged 2.9 
percent of the mean diameters for both species and all size classes. The o for needles 
was obtained from measurements on needle cross sections (Brown 1970). 


Depth.--Twenty-four measurements of depth were systematically taken just prior to 
ignition of each plot. Standard errors of the estimate ranged from 1 to 5 percent of 
the mean depths. Depth values were used to compute packing ratios. The packing ratio 
is the ratio of fuel volume-to-volume occupied by fuel and it indicates compactness of 
the fuel for internal calculations in the spread model. 


Ait 
hj 


ial NAY 


Figure 2.--Variations existed tn sizes and arrangements of particles, as shown for 
ponderosa pine (top) and Douglas-fir (bottom), even though all plots were loaded 
as untformly as posstble. Clumps of particles extst because of the nature of 
branching and foliation. 


Particle density.--Density values were determined from 12 randomly picked samples 
from each species and diameter size class. A mercury pycnometer furnished the density 
values for the needles and 0- to l-cm. size classes. For the other size classes, 
volumes for calculating density were determined using average diameters and a length of 
15 cm. for the samples. Average diameters were based on two measurements taken perpen- 
dicularly at l-cm. intervals. Weights were on an ovendry basis and volume on an air-dry 
basis. Standard errors of the estimate were less than 4 percent of the mean density 
values for each size class. 


Low heats of combustton, total ash content, and stlica-free ash content.--Duplicate 
determinations for three samples picked randomly from the slash provided these values 
for needles and the O- to 1- and 1- to 5-cm. size classes. 


Motsture content.--Six samples for each size class, three near the top and three 
near the bottom of the slash, were collected just prior to ignition. The branchwood 
was ovendried and the needles subjected to xylene distillation. Standard errors for 
each plot ranged from 3-1/2 to 5 percent of the mean moisture contents. 


Wind.--Wind velocity was measured at the expected midflame height close to each 
plot but away from indraft influences of the fires. Anemometers (having a threshold of 
0.5 m.p.h.) connected to portable battery-operated recorders furnished the measurements. 
A wind vane also placed at midflame height recorded wind direction. 


Temperature and relative hwnidity.--Temperature and relative humidity were measured 


just prior to ignition and continuously. They were recorded for reference but were not 
required for calculations in the spread rate model. 


Measurement of Fire Variables 
The plots were ignited across one end and allowed to burn for 5 feet to attain a 
quasi-steady state before spread rate measurements were begun. Time for the fire to 


burn between 2-foot intervals was recorded over a 20-foot distance (fig. 3). 


Reaction intensity was determined from: 


a (c n,w hs )/t (2) 
i=l 
T= DAR (3) 
where: 
tT = reaction time of flame front (time for propagating flame to pass a fixed 


point), min. 
n = average fractional weight loss of fuel in flame front 
w = loading, 1b./ft.? 
h = low heat yield, B.t.u./1b. 
i = index for particle size classes 
D = flame depth, ft. 


R = rate of fire spread, ft./min. 
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fire to pass each of the markers shown, which were spaced 2 feet apart. 


Ovendry weights, which were obtained before and after the burns, permitted calcula- 
tion of the particle weight loss within the propagating flame front. Weight loss by 
particle size class (n) was estimated at one location in each plot. Two 15-cm.-long 
sample branches per size class were tagged and located in a 20- by 100-cm. area 1.5 m. 
from the end of each plot. The fires were quickly extinguished using a water fog after 
the rear of the propagating flame passed over the particles. 


Flame length and depth were measured on photographs of the fires' profiles taken 
at 2-foot intervals along the plots. Flame length is the distance. from the outermost 
tip of the unbroken flame to the center of the flame base which lies at the top of the 
fuel. Flame depth is the horizontal distance through the base of the flame from the 
leading edge of the flame to a point at the rear where it tapers rapidly into a low 
and often sporadically flaming region. 
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RESULTS AND DISCUSSION 


Good control over the fuel variables (see tables 2 and 3) was maintained by care- 
ful construction of the plots and by obtaining high precision in estimating them. 


Rate of Fire Spread 


The mathematical model predicted rates of spread that were higher than observed 
rates on all but three of the plots (fig. 4 and table 4). Deviations (predicted values 
minus observed values divided by observed values) ranged from -14 to 580 percent. The 
average overprediction (percent deviation) was 82 for the ponderosa pine plots and 121 
for the Douglas-fir plots. These are in reasonably close agreement considering the 
complexity of the mathematical model and the slash fuels. 


Several reasons beyond experimental error probably account for the deviations. 
The consistent overestimate of the spread rates indicates bias that could be explained 
by certain aspects of the mathematical model: (1) Assumptions about the fuel; (2) the 
moisture of extinction chosen for the slash; and (3) the method of weighting input 
parameters by fuel surface area. 


Some bias may have been introduced on the plots having higher intensity fires 
(some of the 30-ton-per-acre plots) because the 8-foot plot width limited flame develop- 
ment. However, this bias should be small. An analysis of the effects of plot width 
on flame size and related rate of spread indicates that an 8-foot plot width would re- 
tard only slightly the rate of spread for fires having flames greater than 4 feet in 
length (Anderson 1968). No retardation of spread rate was indicated for fires having 
flames less than 5 feet in length. 


ata 


Table 2.--Average fuel properties of the slash according to spectes and stze class 


Fuel : Ponderosa pine : Douglas-fir 
roperty > Needles i 0s a asi 0-1 13 ses 


Low heat valuel/ 
(Bt. / ib.) 8,826 9,274 8,885 8,885 8,898 8,691 8,691 


Silica-free ashl/ 
(fraction of dry 
weight) .0155 .0123 . 0086 .0086 .0145 . 0097 .0097 


Total ashl/ 
(fraction of dry 


weight) .0387 .0155 0096 0096 .0166 0108 0108 
Particle density o/ / 

(1b./ft. 3) aes 35.6 33.0 30.8 ELSA See SAEs 3205 
°(£t.2/£t. 3) 1,756 156 73.5 32.9 393 80.7 29.6 
Proportion of total Sy) 3/ 

loading (percent) 22 i 43 28 = 27 = NS 347 


1Branchwood for the 1-3 and 3-5 cm. size classes was combined for the chemical 
analyses. | 

2Density values for the 0-1 and 1-3 cm. size classes were not significantly differ- 
ent using a Tukey test at the 0.05 confidence level. 

3These values were not used to determine loadings for the spread calculations. 
Instead, loadings by size class were sampled for each plot. 
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Naas line 


Legend 


© Ponderosa Pine (15 T/A) 
@ Ponderosa Pine (30 T/A) 
O Douglas Fir (15 T/A) 
™ Douglas Fir (30 T/A) 


Predicted rate of spread (ft./min.) 
o 


Figure 4.--Comparison 
between observed and 
predicted rates of 


spread for all plots. 0 2 4 6 8 10 12 
Observed rate of spread (Ft./min.) 


*pouTquod 919M TOADT BDUEpTzUOd 
S0°O 942 3B 380} adUeTeFFTP JUBSTFIUBTS ATIsouoY s,AayNJ BuTsn QustazFtp ATIUVITFTUBTS JOU O1OM JeYI SSETD OZTS Yee LOZ SUBeW, 


88 

T1z 
8ST 
S8T 
Trt 
Z£T 
Laat 
vit 
022 
88 

SLb 
I1z 
TLC 


PST 
46 
9LT 
901 
80¢ 
LOT 
L6 
90T 
OL 
OL 
S8T 
OST 
TIz 


Un) * 44 


poodsputm 


vz 18 009T 6S" 
8e 3/6 OOTT So'T 
))| 18 0002 TS 
82 9L 0002 vol 
99 Ss 0£ 80 82'T 
S8 bp 0020 oo 1 
1Z v8 ST6T 76° 
LZ SZ S102 62a 
4 SZ 0002 67° T 
bz 48 006I 9ST 
L£@ 69 OL9oT $9° 
LZ 78 OL6T Ov'T 
0¢ gL OTTT (Whe 
49 oP 0080 14°ty 
LL 1S 0080 59° 
LZ SZ 00TZ £6° 
s¢ 82 0002 Jeni VI 
9¢ $2 OOTT YAY I 
09 £S 2:10) 78° T 
£8 9b 0040 76°Z 
8¢ bl 0002 $9” 
0S 89 O00TZ L40°T 
8¢ s9 O00TZ LO°T 
bl 8b 0080 L9eT 
9L zs S¢2o gio 
(4) 9S 2:10) yl°2 
quaodad "do "aH "hd 
AqYtptumy ainjzesrodua y Out y, yadop 
OATIBTOY pony 


mmamouwUuvoCooanntan 
ADDON NMDOAOMAAS 
[on io! ei 


a els or 
ce 


. 


or~stsMNononrnnstoo 
omonnnonnoanna- 


Lo | 


. 
. 


et 


i- 


. 


AMANDOUOOMUOAADA 
ADWAAnMNNHAOHRAGS 


AmANUDUOUOTODOANDA 
AAAANMmODODAONUOAD 


Loni 


. 


Oon~neryrTNONwemarNNAOM 
omonm~nmnmonrnOoOMnNaeo 


o~neruNmMowmaonnnwaon 
onm~aonwouowueonnmonane 


Co~mOMNMOoMnNnNNnNANS 
ernmronMnvo onnnuonve- 


* €-T * 1-0 ‘ aTpeeN 
I ue {U0D 91n stow 


£8" 9c” Ob" Ey 
SOR oe 09° [be 
vs" Ge cb (A) 
VOR, OG 62° Sia 
LO L9G 99% ve* 
SOME SCL: £9" Loe 
s8° 1 bb [ee 
£8° 6£° 9os 8° 
£8° 6g" 9 8t° 
SOT OOM ies 2S og * 
z8° Sf vA 365 
SOT SOs (ZS Tes 
18° (ts gt" Lee 


Tel ve ES 60° LG 
bs° vz" oo 90° 8I° 
ov TI ih: A} Tite (AX 
Las 0z* o£* S0° Sis 
Sie Oa [9 oe Eki 
La 02° 0° s0° oT” 
ov Ten OPs [95 10) (hie 
Tay Okay 0s* s0° oT 
Spel “OV” Ss 0) (4 
We 07° (0} 9% so" OL” 
lve 0%. 19° OT’ (A 
12° 02° 10 S0° ois 


CP TORO T9eOm ORO" Faso 


TBIOL' S-¢ * S-T * T-0 * @eTPeeEN 
But peoy 


Buqurng fo auny ay, 4D 407d yova uof suo1gzpuoo 7HZUaUUOU2LAUAa puv saiquadoud 7anj--"*¢ 9TQeL 


EC 
14 
eZ 
02 
81 
et 
oT 
qst 
esl 
bl 
£T 
4) 
c 


92 
£7 
Iz 
Un 
OT 


ATF 
-se,3noqg 


outd 
eBsoLopuog 


satoeds 


13 


Table 4.--Rate of spread, flame length, and flame depth for 
tndtvidual plots 


Plot Rate of spread Flame ey 
Observed ; Predicted donee crepe 
Ft. /min. Ft. /min. Mice Er 


PONDEROSA PINE 


1 ORS 1525 Bite) Se 
2 1S 4.71 Ib) ew 
3 Sais 5.00 5.4 4.2 
4 a e)il Weis: ZZ Wed 
6 Oy 1.48 7.4 BS) 
7 .61 3.0 Weal 0) 
8 4.14 (Ose IES) 6.1 S58) 
9 1.86 6.34 WoW) 1.4 
10 Posy 9.47 a5) Joo 
ill .90 ie (/0) cobs) Boe 
21 eS ZS Sho Boo) 
23 .78 Om 2.0 eS) 
26 1.54 INt3)5) S55 Zine 
DOUGLAS-FIR 

5 . 34 Zo a) 4 
2 2.63 Swe) 4.9 Se 
NS 1.68 2.34 aI) of) 
14 1.68 Dishes 4.0 Sia 
15a Zd9 3.86 2.4 Lets) 
15b 1s ZO Was) Tye) 
16 55 2.84 ied og) 
Ie .87 PHAR) 3.4 2.4 
18 84 Deas) M5 (0) 0 
20 2.88 2.48 Toe 4.4 
22 lez Wes 2.6 Le) 
24 La 52 2.43 By) 1.4 
25 0 Ili 0 0 


1Depths were determined from measurements on photographs except 
for plot 24 which was determined by a visual estimate. 
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Figure 5.--Sporadic flaming typtfted the 15-ton-per-acre plots; this was probably caused 
by the relatively small amount of fuel avatlable and the clumpiness of particles. 


Discontinutty of fuel and heat flux.--Deviations were highest for the lowest values 
of packing ratio and fuel loading (table 5). This was especially true for the ponderosa 
pine plots. A discontinuity in the arrangement of fuel particles and a corresponding 
decline in the supply of heat for preignition could have existed because of the spac- 
ing between fuel particles. The spacing was greater in the least compact fuels: the 
branching habit and manner of needle growth created various sized gaps between particles, 
even though the slash was distributed uniformly during construction of the plots (fig. 
2). Some gaps were probably so large that a substantial propagating heat flux was re- 
quired to ignite unburned particles. In contrast, the mathematical spread model was 
developed primarily using fuel beds having evenly spaced particles without any large 
gaps between them. 


Deviations for the plots having the highest loadings were less, possibly because 
the higher fire intensities furnished ample heat flux for a nearly uniform rate of 
particle ignitions. Fires in the low loading plots generally burned more sporadically 
because particle ignitions were not occurring at a uniform rate (fig. 5). 


On one of the low loading Douglas-fir plots (no needles present) fire would not 
sustain itself; on another similar plot, fire spread was barely sustained. The amount 
of fuel, particle size, and particle spacing were the dominant limiting factors. The 
spread model contains mathematically continuous functions that are defined by scant 
data at the limiting end of the fire scale. 
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Table 5.--Percent devtations for the observed and predicted spread rates 


Loading 
tons/acre 


15 
15 


15 
aS) 


TS 


15 
Average 


25 


Average 


30 
30 


30 
30 


30 
30 
Average 


Packing ratio 


0.014 
012 


.021 
.019 


-034 
.040 


.016 
.015 


.027 
.026 


.041 
.049 
.043 


.019 
.019 
O27, 


.038 
.032 


.041 
. 048 
.035 
031 


.040 
.039 


.047 
.047 


Percent 


PONDEROSA PINE 


DOUGLAS-FIR 


169 
241 


117 
200 


113 
30 
145 


al 
48 


33 
25 


11 

25 
] 

34 


segregated by spectes, loading, and packing ratto values 


Deviation 


Average percent 


205 


158 


72 


60 


29 


12 


95 


206 


22 


54 


178 


37 


lomitted because no observed rate of spread resulted in an infinite 


percent deviation. 


16 


Moisture of extinction.--Rate of spread is directly proportional to the moisture 
damping coefficient (n,), which enters the mathematical model as a multiplier of reac- 
tion intensity. 


a= a= 2 is9 pits ee 28 Lig S59 ue (4) 
; M : M i M 
x x x 
where: 
Me = moisture content of fuel, fraction of dry weight 
M. = moisture of extinction, fraction of dry weight 


The moisture of extinction is the moisture content of fuel at which fire cannot 
sustain itself. Criteria for choosing M, are not well established; thus use of this 
function in the model is rather subjective. We used an My, of 0.24. Working backwards 
through the model using fuel inputs and measured rates of spread for each plot, we 
found that an average M, of 0.11 for the ponderosa pine and 0.13 for Douglas-fir would 
have provided perfect agreement between observed and predicted spread rates. My, values 
that permitted perfect agreement ranged from 0.065 to 0.45 for all plots. There is good 
reason to suspect that an M, of 0.24 is too high for the fuels tested. If so, this 
would account for much of the deviations between observed and predicted spread rates. 


A comparison of M,. with packing ratio and loading values suggests that a relation- 
ship exists between My and these fuel properties. In this study, plots with higher 
loadings developed more intense fires. The M, would be expected to increase at higher 
intensities because more energy is available for preheating fuel to ignition. Plots 
having lower packing ratios (more porous) had greater percent deviations between 
observed and predicted spread rates. Apparently this was in part caused by discontin- 
uities in fuel and heat flux. Discontinuities should tend to lower M,. 


Surface area influence.--Weighting the input variables, especially o, by the amount 
of surface area in each particle size class possibly leads to overprediction of spread 
rate. The thinnest fuel particles in a fuel complex receive heavy weighting for their 
influence on fire spread rate. This seems appropriate, but perhaps too much weighting 
is received. Further study should refine the accuracy of this function in the model. 


Regression Analysis for Fire Spread 


A step wise multiple regression analysis using the combined data for both species 
of slash resulted in the following equation: 


Y = 0.109(X,) + 2.01(X2) + 0.00827(X3) - 1.96(X,) (5) 
where: 


Y = rate of spread, ft./min. 


X1= surface area of fuel 0 to 1 cm. in diameter per square foot of ground, 
dimensionless 


Xo= total fuel loading, 1b./ft.? 
X3= windspeed ft./min. 

X,= bulk density of fuel, 1b./ft.? 
R2= 0.71 


Standard error of estimate (0.05 confidence level) = 1.04 


Zs 


Surface area of the O- to l-cm. fuel accounted for differences between species. 
Presence of needles on the ponderosa pine branches provided this species with the most 
surface area. Of the correlation between the X and Y variables, X, was closest to 
Y (r = 0.57). Moisture content was nonsignificant in the regression because of the 
small range in moisture values. Xj was nonsignificant in regressions in which species 
were kept separate. 


Reaction Intensity 


The measured IR values were mostly less than the predicted ones, as shown in the 
following tabulation: 


Ponderosa pine Douglas-fir Speetes combined 
Average reaction intensity 
Measured, B.t.u./ft.?-min. 4,621 3,599 4,130 
Predicted, B.t.u./ft.?*-min. 4,800 3,808 4,324 
Average absolute deviation, percentl/ 53.6 52.8 Sore 
Overpredictions 10 9 19 
Underpredictions 5 4 7 


We could not substantiate conclusions about the accuracy of the model to predict Ip 
because we were unable to measure Ip precisely under field conditions. There was a 
tendency to overestimate flame depth used to calculate Ip (equation 3) because the 
widest portion of the flame was outlined on the photographs. The fact that fuel weight 
loss during passage of the flame front was based on an average value for all plots also 
contributed to the deviations. 


Weight Loss Related to Particle Size 


The average percent weight loss during passage of the flame front for sample fuel 
particles from all plots is tabulated below: 


Ponderosa pine Douglas-fir 
Percent Percent 
Needles 2/ 96 -- 
Branches by) 
diameter’ (cm. ) 
0.2 = 4 oq 4 2 
Oz5 -- 80 + 16 
1 (RH as ANS Ws cewalfs) 
1-3 A202) 364 15 
3-5 iss es 1/ M6 


IMeasured Ip - predicted Ip x 100 was calculated for each plot and averaged with 


Measured Ip 

sign ignored. 

The needles were completely consumed; however, needles contain 4 percent ash which 
theoretically is not a weight loss. 

3Diameters were chosen to provide sensitivity at the small sizes. Inasmuch as 
ponderosa pine and Douglas-fir branches have a different size distribution, different 
diameters for the small sizes were chosen. 

"The second values are standard deviations. For ponderosa pine particles, n = 26; 
for Douglas-fir particles, n = 18. 
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The percent of total fuel weight loss accounted for by needles and branchwood in 
the different particle size classes is shown in the following tabulation:9/ 


Ponderosa _ pine Douglas-fir 


Pereent Percent 
Needles 44 -- 
Branchwood by 
diameters (cm. ) 

0-0.5 -- 19 

0.5-1 a1 25 

1-3 Si7/ 43 

3-5 8 13 


Percentages for individual plots did not vary significantly for different loadings or 
bulk densities. 


When the above percentages were determined, we knew the percent weight loss of 
individual particles and the proportion of the total fuel contributed by the various 
size classes. Fuel less than 1 cm. in diameter accounted for half of the total weight 
loss even though this size material was less than 30 percent of the total loading 
(table 1). Fuel greater than 3 cm. in diameter provided only about 10 percent of the 
total weight loss in the flame front (fig. 6). No change in diameter was observed for 
pavercles of this, ‘size. 


For slash and other fuels containing a similar mixture of particle sizes, our 
findings indicate that only the fine fuel components, essentially particles less than 
3 cm. (about 1 inch in diameter) supply the energy that characterizes propagation of 
the spreading flame front. A generalization of this statement for all fuels likely 
would be incorrect. Different proportions of fine fuels probably correspond to dif- 
ferent diameter limits for the component contributing most of the energy to the propa- 
gating flame. 


The percent weight loss of particles is compared in figure 7 with the effective 
heating number of the mathematical model defined in equation 1.&/ The effective heating 
number is a function of particle size; it expresses the percent of a particle that is 
heated to ignition ahead of the fire. In figure 7, the percent weight loss curve is 
drawn through points that represent particle size and the average percent weight loss 
in the propagating flame front for the sample particles. 


The relation between the two curves, each representing a different variable, was 
expected and supports the concept expressed by the effective heating number. The two 
curves are close together at the very small particle sizes because the small particles 
are heated throughout at ignition and their entire organic mass rapidly converted to 
heat energy. As particle size increases, up to some point, the amount of a particle 
required for heating to ignition (e€) becomes less than the amount converted to heat in 
the flame (n). The effective heating number only involves that part of the particle 
that receives heat during the time when the particle surface rises to ignition tempera- 
ture. Weight loss in the flame front, which takes place primarily after ignition, is 
characterized by rapid heat transfer and combustion, and thus involves a larger propor- 
tion of the particle. The percent weight loss (on an ash-free basis) should always 
exceed the effective heating number for the same particle sizes. | 


“The percentages are averages based on all plots. 
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Figure 6.--After the flame front passed here, the fire was extinguished with a fog 
nozzle. Most of the remaining branchwood was over 1 em. tn diameter. About 80 
percent of the branchwood less than 1 em. tn dtameter was transformed into heat 


energy. 
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Legend: 

8 80 @ Ponderosa Pine 80 
2 @ Douglas Fir 5 
= x Effective heating number ‘= 
= 60 60") 0 
re) mo 
ie & 
wre ce] 
D O 
in AO AOais 
> = 
o 
i= 2 
S a 

2510) 20 

0 

@) .005 .010 O15 020 025 .030 


i/o - ft 


Figure 7.--Percent weight loss in the flame front and percent heated to ignition as a 
funetton of the particle surface area-to-volume ratto. 
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Flame Length 


The relationship between flame length and certain fuel properties was examined 
using a stepwise multiple regression analysis. The fuel properties that varied sub- 
stantially from plot to plot served as the independent variables. The following equa- 
tion resulted using data from all plots: 


Y = 4.28(X1) + 0.00151(X5) - 2.92 


Y = flame length, ft. 


Sj) — tatal loading, 1b./ft.2 

Xp = average particle surface area-to-volume ratio weighted by fuel surface 
area in each size class, ft,-1 

R* = 0.60 


Standard error of estimate (0.05 confidence level) = 1.52 


The packing ratio was nonsignificant as an independent variable. Flame length was 
expected to correlate with packing ratio because (1) flame length relates closely to 
rate of energy release (Rothermel and Anderson 1966; Thomas 1963), and (2) packing 
ratio is quite sensitive to rate of energy release for the smaller size fuels in the 
spread model. If data covering a greater range of packing ratios had been acquired, it 
probably would have been significant in the regression analysis. However, loading of 
fine fuel (certainly particles less than 2 to 3 cm. in diameter) appears to be by far 
the most important physical property relating to flame length. 


6W. H. Frandsen. The effective heating of fuel particles ahead of a spreading fire. 
USDA Forest Serv., Intermountain Forest and Range Exp. Station, Ogden, Utah. (Manuscript 
in preparation. ) 
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CONCLUSIONS 


The mathematical model's predictions correlated reasonably close to observed rates 
of spread, although the model tended to overestimate. The study showed that M,, dis- 
continuities, and possibly the method of weighting by surface area in the mathematical 
model need additional study. The My appears to be a function of particle size, loading, 
and packing ratio. Its functional relationships in living fuels probably involve other 
properties as well. Weighting the input parameters by surface area in each size class 
seems appropriate for heterogeneous fuels. However, the contribution of the very fine 
fuels to spread rate may be overemphasized. 


Although imprecise field methods prevented extensive analysis of I, and related 
flame length, refinement of Rothermel's model should improve their prediction. Con- 
tinued development of mathematical models for fire spread and intensity should be pur- 
sued because they will provide a sounder basis for quantitative fuel appraisal. 
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